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ABSTRACT 

Stars in disks of spiral galaxies are usually assumed to remain roughly at their birth radii. This 
assumption is built into decades of modelling of the evolution of stellar populations in our own Galaxy 
and in external systems. We present results from self-consistent high-resolution A^-body -I- Smooth 
Particle Hydrodynamics simulations of disk formation, in which stars migrate across significant galac- 
tocentric distances due to resonant scattering with transient spiral arms, while preserving their circular 
orbits. We investigate the implications of such migrations for observed stellar populations. Radial 
migration provides an explanation for the observed flatness and spread in the age-metallicity relation 
and the relative lack of metal poor stars in the solar neighborhood. The presence of radial migra- 
tion also prompts rethinking of interpretations of extra-galactic stellar population data, especially for 
determinations of star formation histories. 

Subject headings: galaxies: evolution — galaxies: spiral — galaxies: stellar content — Galaxy: solar 
neighborhood — Galaxy: stellar content — stellar dynamics 



1. INTRODUCTION 

In a universe where mass assembles by accretion of pro- 
gressively larger constituents, thin disks of galaxies form 
during quiescent evolution following the major m erging 
epoch (jBrook et all 120041 : [Robertson et al.l l2006f ) . The 
inner parts, with lower angular momentum assemble 
first, followed by the higher angular momentum compo- 
nents, resulting in "inside-out" build up of disk material 
(lLarsonlll976HWhite k, Frenklll99lHMuiioz-Mateos et all 
l2007t ). Because gas can efficiently dissipate energy, its 
motion around the galaxy is mostly circular. Stars born 
in such gas disks are initially on nearly circular orbits, 
but their infant kinematic state is highly fragile; non- 
axisymetric perturbations such as bars or spiral arms 
readily drive the orbits away from circular. Despite the 
increase in eccentricity, mean orbital radii are assumed to 
remain relatively constant. As a result, interpretations 
of Galactic and extra-galactic stellar population observa- 
tions invariably make the fundamental assumption that 
stars remain at roughly their birth radii throughout their 
lives. Consequently, this assumption is built into the vast 
majority of models of formation and evolution of galac- 
tic disks over the pas t few decad es (e.g. ITinslevI 119751: 
Matteucci fc FrancoislligsllCarigill lOQetlChiap pini et all 



1997t iBoissier fc Prantzog|1999[ ). 



Recent theoretical and observational evidence chal- 
lenges this static picture. Spirals have long been known 
to be an important source of kinematic heating in galac- 
tic disk s, gradually increasing th e eccentricities of stellar 
orbits (jJenkins fc BinnevI 119901 ). However, if a star is 
caught in the corotation resonance of a transient spiral 
it may move inward or outwar d in radius while preserv- 
ing the circularity of its orbit (jSellwood fc Binnevll2002l 
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hereafter SB02). Transient spirals with a wide range of 
pattern speeds are present throughout the evolution of 
the disk. Different pattern speeds result in spatially dis- 
tinct locations of corotation resonances, allowing a star 
to undergo a series of resonant encounters in its lifetime, 
riding the spiral waves across large portions of the galaxy 
while retaining a mostly circular orbit. Stars born in- 
situ may therefore remain kinematically indistinguish- 
able from those that hav e come across the galaxy. 

In lRoskar et al.l (|2008f ) we presented first results from 
our A^-body Smooth Particle Hydrodynamics (SPH) 
simulations of disk galaxy formation in which stars 
migrated radially due to their resonant interactions 
with transient spirals. These migrations yielded ra- 
dial stellar population gradient predictio ns, which have 
recen tly been observed in resolve d-star (Ide Jo ng et ail 
2007f). and integr ated- light studies (jAzzoUini et al. 2005 



Bakos et al.l2008[ ). In conjunction with our models, these 



observational data strongly imply that radial migration 
is an important process in observed galactic disks. In 
this Letter, we present fur t her an alysis of the simulation 
presented in iRoskar et al.l (|2008D . focusing on the wide- 
ranging implications of radial migrations for the observ- 
able properties of stellar populations. 

2. SIMULATION 

The initial condi tions are created as in 
iKaufmann et al.l ()2006l ). and are designed to mimic 
the quiescent stage following a last major merger when 
thin disk formation commences. The s ystem consists 
of a spherical dark matter NFW halo (jNavarro et al.l 
Il997t ) in which we embed a spherical halo of gas with 
the same initial profile. The gas halo is initially in 
hydrostatic equilibrium. The total mass of the system is 
IO^^Mq, analogous to a Milky Way-type spiral galaxy, 
with the baryons contributing 10% of the mass. We 
resolve the system with 10^ particles per component, 
resulting in initial mass resolutions of lO^M© and 
IO^Mq for dark matter and gas respectively. Star 
particles form with masses that are a fraction of the gas 
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particle mass, resulting in typical star masses around 
SxW^Mq. To form a rotationally-supported disk, we 
also impart angular momentum to the gas component 
corresponding to a spin parameter value A=0.039, 
as r notivated by cosrn ological iV-body experiments 
(e.g. iBuUock et al.l l200lh . We evolve the system using 
the A^-bod y + Smooth Partic le Hydrodynamics code 
GASOLINE IWadslev et al.ll2004D for 10 Gyr. 

Once the simulation begins, the gas cools and collapses 
to the center of the halo, forming a thin rotating disk 
from the inside-out. When the gas reaches densities and 
temperatures which allow for star formation, the star for- 
mation and superno va (SN) feedback cycles are initiated 
(jStinson et al.ll2006D . Since our disks form without any a- 
priori assumptions about the interstellar medium (ISM) 
or the stellar populations present in the disk, we can 
self-consistently follow the evolution of their properties 
as the disk grows. Transient spirals cause radial redistri- 
bution of stellar material in a manner analogous to the 
mechanism presented in SB02. Due to space limitations, 
a detailed discussion of these processes is deferred until 
a future paper. 

Our simulations do not account for the effects of evo- 
lution in a full cosmological context, but in the standard 
ACDM paradigm mergers are significantly more impor- 
tant at early epochs. The simulation presented here is 
therefore designed to mimic the quiescent phase of disk 
galaxy formation during which the thin disk forms. The 
increased resolution gained by stepping outside of the 
cosmological context allows us to isolate the effects of 
key dynamical processes, which determine multiple ob- 
servational characteristics. Despite these simplifications, 
physical properties of our model such as the rotation 
curve, star formation rate, disk scale-length, and disk 
mass fraction, are analogous to those of observed sys- 
tems. Pre-enrichment of merger fragments would not 
change the details of disk build-up, but simply offset the 
metallicity distribution. 

3. RADIAL MIGRATION AND ITS IMPLICATIONS 

Models of galactic chemical evolution have been enor- 
mously successful in explai ning the properties of stars 
in our solar neighborhood (iMatteucci fc Francoiill989t 
Carid 119961 : IChiappini et al.lll997t iBoissier fc PrantzosI 



19991 ). One example of such model results is the age- 
metallicity relationship (AMR). The AMR is expected to 
arise due to progressive enrichment of the ISM through 
stellar feedback. Stars of the same age in the same 
general region of the galaxy are therefore expected to 
have similar metallicities. Indeed, early determinations 
of the AMR confirmed that the mean trend of stars in the 
solar neighb orhood is towa rd lower metallicity with in- 
creasing age dT waro i figsl . Models, which assume that 
stars remain where they are born and return their nu- 
cleosynthetic yields to their local ISM, typically success- 
fully reproduce this trend. However, evidence suggests 
that a large amo unt of scatter is present in the AMR 
of fie ld stars fjEdvardsson et al.l 119931; iNordstrom et al.l 
|2004[ ) and open clusters (jFriel et al.ll2002l ). meaning that 
stars of the same age are observed to have a wide range of 
metallicities. In the framework of a galactic disk where 
radial annuli arc self-enriched, homogeneous, and iso- 
lated from one another, old metal-rich stars or very young 
metal-poor stars are an impossibility. 



Allowing stars from slightly different regions of the 
disk to enter the local sample due to the eccentricity of 
stellar orbits has been considered as a possible explana- 
tion, bu t can account only for up to 50 % of the observed 
scatter (Nordstrom et al. 2004; Binnevl l2007l) . The large 
amount of scatter in the AMR therefore poses an impor- 
tant challenge for models of disk formation and should 
be considered as large a constraint as the mean trend 
(jCarraro et al.i ri998). The unexpected degree of scatter 
implies that stars in the solar neighborhood were either 
formed from a highly inhomogeneous ISM or have come 
from wildly different parts of the galaxy. 

In our simulation, the latter option offers an enticing 
solution. In the left-most panel of Figured] we show the 
distribution of birth radii for stars, which at the end of 
the simulation are on nearly circular orbits within an 
annulus analogous to the solar radius. We define 'so- 
lar radius' as a general region of the disk interior to the 
disk break, but approximately 2-3 scalelengths from the 
center. In our simulated galaxy this region is between 
7-9 kpc (indicated by dashed lines in Figure [T]). The 
black line represents all stars while the blue and red lines 
show the distributions of metal poor ([Fe/H] < —0.3) and 
metal rich ([Fe/H] > —0.1) stars respectively. Roughly 
50% of all "solar neighborhood" stars have come from 
elsewhere, primarily from the disk interior. Interestingly, 
some metal poor stars have been scattered into the so- 
lar neighborhood from the outer part of the disk. Such 
migra tion has recently been inferred from observational 
data (jHaywoodI l2008 (i. Metal-rich stars, like our Sun, 
could have originated almost anywhere in the Galaxy. 

3.1. Age-Metallicity Relationship 

It has previously been suggested that such large ra- 
dial migrations could account for the scatter in the ob- 
served AMR (SB02, Binncy 2007|), but our models al- 
low us to examine the implications of migrations in a 
self-consistent simulation of a growing and star-forming 
galactic disk. We show the AMR as determined from our 
simulation in the middle panel of Figure [TJ Colors indi- 
cate relative particle densities for clarity in high-density 
regions. Diamonds and error bars indicate the means 
and standard deviations in several broad age bins. With 
the exception of the oldest age bin, the AMR is fairly 
flat. The relatively high number of very old (> 9 Gyr) 
extremely metal-poor stars is partly a remnant of our 
initial configuration, which includes only pristine, metal- 
free gas. This is likely to be different in a cosmological 
setting where a fraction of the gas is accreted from sub- 
haloes that have already formed stars and been enriched 
by stellar feedback. Also, note that we assume each star 
is born with the mean metallicity of the ISM at its birth 
radius so that we can more easily isolate the effects of 
radial migration on the various observational relations. 

To illustrate the impact of the dynamical evolution of 
the disk, we plot the AMR in these same age bins by us- 
ing only stars that actually formed in the solar neighbor- 
hood. This latter relation is analogous to that derived by 
models which assume a dynamically non-evolving disk, 
essentially only depending on the nucleosynthetic yields 
and star formation and gas infall rates at the solar cir- 
cle. In comparison to the full simulation, the in-situ re- 
lation shows a much tighter correlation between metal- 
licity and age. Inclusion of full dynamical modeling sig- 
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Fig. 1. — Properties of stars in the solar neighborhood. Left: Histogram of birth radii for stars that end up in the solar neighborhood 
on nearly circular orbits. The black, red, and blue lines represent all, metal-rich, and metal-poor stars respectively. Middle: The age- 
metallicity relation: color contours represent relative particle densities where point density is high. Diamonds and error bars indicate mean 
values and dispersion respectively. Squares show the AMR if stars are assumed to remain in-situ. A small horizontal offset is applied to the 
two sets of symbols for clarity. Right : Metallicity distribution function (MDF): t he si mulated distribut i on is shown with the solid black 
histogram; diamonds and asterisks show data from IRocha-Pinto fc Maciell I I1996I) and IHolmberg et al] II2007I ) respectively. The dashed 
histogram is the MDF if stars are assumed to remain in-situ. 

nificantly flattens the AMR and increases its dispersion 
by at least a factor of two. We note that although our 
simulations are not designed to reproduce Milky Way 
observables, they successfully reproduce the major qual- 
itative features of the observations: the observed AMR 
is indeed relatively flat with a high degree of scatter ^ 



stars 



Gas 



IS maeea relativel y iiat witn a nigii deg ree oi scatter 
(jEdvardsson et all Il993l : iNordstr om et al.l |2004) . which 
gradually increases with age fHavwoodl I2008D. This 



gradually increases with age 
agreement strongly suggests that a large amount of radial 
migration is the missing piece of the AMR puzzle. 

3.2. Metallicity Distribution Function 

If radial migration is important for the AMR, then it 
should also leave a significant imprint on the metallicity 
distribution function (MDF), another key observational 
constraint for models of Galactic chemical evolution. The 
simple, closed-box chemical evolution picture where no 
material is allowed to enter or leave a given radial annulus 
has long bee n recognized as inade quate at explaining the 
local MDF (|van den Berghl ll962l. However, the dearth 
in the relative number of low- metallicity stars, known as 
the G-dwarf problem, is usually explai ned b y allowing for 
the in flow of metal-poor gas (,Larson|[l974l : [l7vnden-Belll 
Il975f ). Gas infall "solves" the G-dwarf problem because 
it allows for prolonged star formation, which results in a 
relatively high number of metal-rich stars. 

Our simulated MDF is shown with diamonds in the 
right panel of Figure [1] conipared to ob servational data 
from Rocha-Pinto &: Maciell (119961 ) and IHolmberg et al.l 
(|2007 ). represented by diamonds and asterisks respec- 
tively. We also plot the MDF of stars formed only in 
the solar neighborhood (shown with squares). Compared 
to the observed MDF, the in-situ distribution is much 
narrower. Radial migration introduces metal-poor and 
metal-rich stars into the solar neighborhood from other 
parts of the disk, thereby broadening the MDF. 

Unlike the solar neighborhood AMR, which is highly 
dispersed and therefore difficult to use as a robust guide 
for Galactic chemical evolution models, the MDF repre- 




* Note that to shift the model quantities into the metallicity 
range relevant for the Milky Way, we apply an offset of +0.2 dex 
to the simulated metallicities. We show the observed distributions 
not to claim quantitative agreement, but simply to illustrate the 
ramifications of dynamical effects. 



Fig. 2. — Left: Azimuthally averaged metallicity profile at 10 
Gyr for stars in different age bins (young, intermediate age, and 
old stars shown with solid, dashed, and dotted lines respectively). 
Right : Azimuthally averaged metallicity profile of the gas at dif- 
ferent times in the simulation. 



sents a much more stringent constraint. It is clear from 
Figure [1] that radial migration may significantly alter the 
nature of this fundamental property of the solar neigh- 
borhood. 

3.3. Metallicity Gradients 

A third major constraint on models of disk formation 
and evolution is the radial metallicity gradient. If stars 
are assumed to remain in-situ, then stellar tracers of dif- 
ferent ages should provide information about the ISM 
metallicity from which they spawned. In the left-hand 
panel of Figure Owe show the metallicity gradient as de- 
termined from stellar tracers of different ages. The slope 
of the gradient decreases with increasing age of stellar 
population. In most observational studies, such a trend 
would be interpreted to mean that t he slope of the ISM 
metallicity was shallower in the past (jMaciel et al.ll2005l ). 
In the right-hand panel of Figure H) we show the actual 
metallicity gradient in the ISM at different times in the 
simulation. Rather than being flatter in the past, as im- 
plied by the flat gradient of old stars, the ISM gradient 
was actually steeper and evolved in the opposite sense to 
that indicated by the stellar tracers. Radial migration 
causes more mixing in older populations, creating the 
appearance of flatter gradients at early times and lead- 
ing to a decoupling of stellar properties from their birth 
ISM. 
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Fig. 3. — Star formation history in several broad radial bins. 

3.4. Star Formation History 

Implications of radial migration can be extended be- 
yond our own Galaxy to extragalactic studies. In re- 
cent years, state-of-the-art observations of resolved stel- 
lar populations have begun to probe detailed, spatially- 
resolved star formati on histories in disks of external 
galaxies (e.g. Barker e t al.l 120071: Ide Jong et all 120071 : 
IWilliams et all ,2008.: Gogarten et al.ll2008[ ). ~picallv. 
observed properties of stars in a given region of a disk 
are used to construct a color-magnitude diagram, which 
in turn is used to constrain stellar population synthe- 
sis models to yield a best-fit star formation history and 
there fore provide in formation about the formation of the 
disk (|Dolphinl 120021 ). Again, the underlying assumption 
is that stars observed today are found close to the radius 
of their respective birthplaces. 

In Figure [3] we show the 'migrated' (solid) and 'un- 
migrated' (dashed) star formation history (SFH) in our 
simulation for several broad radial bins. The 'migrated' 
SFH is constructed by considering the distribution of 
ages in the desired region of the disk at the final snap- 
shot of the simulation, much like an observer would do 
given a snapshot of a present day galactic disk. The 
'un-migrated' SFH is the actual SFH, recovered by con- 
sidering the times of formation for all stars formed within 



a given radial bin, regardless of the final radius to which 
they eventually migrated. We show here all stars regard- 
less of the eccentricity of their orbits, but our conclusions 
remain the same even if we restrict our analysis only to 
stars on nearly circular orbits. 

The discrepancy between the 'migrated' and 'un- 
migrated' SFH is large, especially for the outermost re- 
gions of the disk. For example, the 'solar radius' bin (7 - 
9 kpc) could erroneously be deduced to have a flat SFR 
for the past ^ 8 Gyr, when it actually didn't reach cur- 
rent SFR levels until just ~ 5 Gyr ago. Migration must 
therefore be accounted for even in extragalactic studies 
of stellar populations. 



4. CONCLUSIONS 

We have shown that radial migration caused by reso- 
nant interactions of stars with transient spirals has signif- 
icant repercussions on the observed properties of a wide 
range of stellar population systems. The range of impli- 
cations and possible solutions to outstanding problems 
given by radial migration is appealing. Further work is 
required to characterize fully the degrees and timescales 
of migration in spiral galaxies of different physical char- 
acteristics. On the other hand, future observational ef- 
forts such as RAVE, GAIA and LSST wiU provide in- 
tricate data sets that should further constrain the radial 
migration process in the Milky Way. Ongoing and fu- 
ture HST projects such as the ANGST and GHOSTS 
surveys, will require radial migration models in order to 
reconstruct the SFH from resolved-star observations of 
nearby spiral disks and thereby learn about the forma- 
tion history of our galactic neighborhood. 
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